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SAVLOC, COMPUTER PROGRAM FOR AUTOMATIC CONTROL AND 
ANALYSIS OF X-RAY FLUORESCENCE EXPERIMENTS 
by Regis F. Leonard 
Lewis Research Center 

SUMMARY 

A program for a PDP- 15 computer is presented which provides for control and 
analysis of trace element determinations by using X-ray fluorescence. The program 
simultaneously handles data accumulation for one sample and analysis of data from pre- 
vious samples. Data accumulation consists of sample changing, timing, and data stor- 
age. Analysis requires the locating of peaks in X-ray spectra, determination of inten- 
sities of peaks, identification of origins of peaks, and determination of areal density of 
the element responsible for each peak. The program may be run in either a manual 
(supervised) mode or an automatic (unsupervised) mode. 


INTRODUCTION 

The recent development of high-resolution solid-state X-ray detectors has made 
possible the use of X-ray fluorescence as a tool for the detection of trace elements in a 
variety of applications. However, analysis of environmental samples by means of X-ray 
fluorescence generates energy spectra for which quantitative hand analysis is prohibi- 
tively tedious. In addition, the number of samples generated by any comprehensive 
monitoring program requires a speed of data acquisition and analysis which is possible 
only through the use of a computerized system. The present report describes the soft- 
ware designed for control and data analysis of such a system. This report should per- 
mit the reader to operate the system. 


PRINCIPLES OF X-RAY FLUORESCENCE AND CHARACTER OF DATA 

When properly stimulated, the atoms of any element emit X-rays whose energy (or 
wavelength) is strictly dependent upon the atomic number of the emitting atom. Since 



these energies may be quickly and accurately measured, and since the X-ray energies 
of all elements are well known and tabulated, it is possible to identify, with little am- 
biguity, the element responsible for the emission. In addition, the number of X-rays 
emitted is proportional to the number of atoms stimulated, so that after a proper cali- 
bration, quantitative measurements of areal densities are possible. 

The X-ray detectors commonly employed provide data in the form of analog electric 
pulses. Data acquisition consists of the digitizing and storing of these pulses, which 
arrive with random timing at the rate of several thousand per second. A schematic 
representation of the X-ray fluorescence facility is shown in figure 1. A typical X-ray 
spectrum is shown in figure 2. This particular spectrum resulted from excitation of 
a sample by using X-rays from a tungsten X-ray tube. It contains a peak (labelled Rh) 
which resulted from the presence in the primary X-ray beam of a small, known amount 
of rhodium, whose purpose was to provide a monitor of the number of photons incident 
on the sample. 

Data analysis consists of determining from this spectrum the number of X-ray 
counts contained in each peak, identifying the origin of the peak, and converting the 
number to an areal density by using the expression 

N(E) = N.(^\ p x f 2e(E) (1) 

Vd^/p, 


where 


N(E) 

N i 

(dcr/dfi ) E 
p x 

n 

e(E) 


number of X-ray photons of energy E detected 
number of X-ray photons incident on target 

cross section for production of X-rays of energy E by exciting X-rays 

areal density of atoms producing X-rays of energy E 

solid angle subtended by X-ray detector 

efficiency of detector for detection of X-rays of energy E 


In practice, only the quantity N(E) is determined absolutely. The number is as- 
sumed to have the form 


N i “ K Rh N Rh 


( 2 ) 


where Np^ is the number of rhodium (monitor) X-rays detected, so that the areal 
density may be written 


2 



( 3 ) 


3 N(E) 

X K ^M§ E ne(E) 


The overall efficiency function 


F(E) = 



(4) 


is then determined experimentally by using samples of known areal density of a number 
of elements. For elements other than those for which F(E) is measured directly, an 
interpolation is carried out with the assumption that the function varies smoothly with 
energy. The areal density then, in terms of directly measurable quantities, becomes 

p =N(E)*® (5) 

N Rh 

Analysis then requires only the determination of the intensity and energy of each line in 
the X-ray spectrum and the intensity of the line used for normalization Np^- 


DATA ACQUISITION AND ANALYSIS 

The entire X-ray fluorescence operation (sample changing, data acquisition, and 
data analysis) is carried out under the control of a PDP-15 computer. Communication 
between the computer and the X-ray fluroescence hardware is by means of a CAMAC 
system (ref. 1). The programmable interface between the computer and the CAMAC 
system has been described in detail elsewhere (ref. 2). The control panel of the 
PDP-15 computer, shown in figure 3, allows the input of a number of parameters as 
well as the selection of a number of different modes of operation, so that a sufficient 
degree of flexibility is provided to make the system useful for a variety of applications. 


Data Acquisition 

The two principal control tasks which the PDP-15 must perform are the changing of 
targets and the control (stop, start, and data transfer) of the 1024- channel analog-to- 
digital converter (ADC), through which all pulse-height data pass. 
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A typical data-taking cycle is initiated upon input of the last piece of sample identi- 
fication data by means of the computer teletype unit. The ADC is turned on, and data 
transfers from the ADC to the computer memory take place autonomously by way of the 
CAMAC interface. During data acquisition the pulse- height spectrum is displayed on a 
video unit, which assures the user that all systems are functioning properly. 

By means of push buttons on the computer control panel, data acquisition may be 
terminated at any time, either by pressing B4, which causes data to be recorded on 
magnetic tape, or by pushing B2, which ends the run without recording of the data. Data 
acquisition also terminates, with all accumulated data being recorded on tape, when the 
elapsed time reaches the preset value as entered (in tenths of seconds) on control panel 
thumb switch TH2. 

Following termination and depending on the mode of operation selected, either 
another run is started automatically (T15 up) or a halt occurs (T15 down) until the oper- 
ator chooses to continue by pressing B2. In either case sample changing may or may 
not be done, depending on the setting of rotary switch R1 on the control panel (sample 
is changed when Rl=2). 

The data recorded upon completion of a run include the run identification (a five- 
character file name and a three-digit file number), the time elapsed during data acquisi- 
tion, and, of course, the pulse-height spectrum itself. 

Because analog- to- digital conversion and data transfers take place at data-channel 
level, they proceed autonomously once initiated, so that the computer central processing 
unit is free at main-stream level to proceed with analysis of a previous data set if so 
instructed. 


Data Analysis 

Analysis of previous data sets may either be initiated by the operator (press Bl) or 
take place automatically when the computer is operating in an automatic mode (T15 up). 
Analysis begins with the recalling of two previous runs from tape: one with the data of 
interest and one that is a background or blank run, which is subtracted to remove events 
resulting from scattering or impurities in the system or sample- supporting matrix. 

When analysis is initiated by an operator, the program requests a file name and 
file number for each of these data sets. In the automatic mode the program searches 
for a data set with the same name as that previously analyzed, increases the file num- 
ber to obtain the next data set, and uses the same background data set. Hence, it is 
always necessary that at least the first set be analyzed manually. 

The first step in the analysis is the determination of the number of events Np^ 
contained in the monitor (rhodium) peak for both the sample data and background data. 
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This is accomplished by locating the rhodium peak and then fitting that portion of the 
spectrum within 20 channels of it with a Gaussian peak plus linear background. The 
number of events within the Gaussian peak is then the of equation (5). The amount 

of background in this region provides the normalization for subtraction of the blank 
spectrum. 

It is assumed that since the blank spectrum results from an unloaded sample of the 
supporting matrix, it is identical in shape to that encountered in the loaded sample. 
Consequently, after normalization and smoothing over three channel segments, the two 
spectra are subtracted channel by channel to arrive at the spectrum to be analyzed. 

The recalled spectrum (with background subtracted) then replaces the live display 
of data currently being accumulated on the video unit. Figure 4 shows the recalled 
spectrum derived from the data of figure 2. The number of events under each peak, 
that is, the number of detected X-ray photons of each energy, is determined by a least- 
squares fitting to each peak of a Gaussian distribution of adjustable height, width, and 
location. Initial estimates of peak locations may be supplied by the researcher by 
means of the light pen and video display (Til down), or they may be located automati- 
cally by the program (Til up). The automatic location procedure is described in ap- 
pendix A. 

Depending on the setting of toggle switches, an energy calibration may be supplied 
at this time (T16 down). The user may supply through the teletype console the energies 
of the first two peaks selected with the light pen. For operation in the automatic mode 
(T16 and T15 up) the program retains the energy calibration from the previous rim. 
Hence, again at least one run must be done manually. The researcher is also asked to 
supply through the teletype console (unless the computer is operating in the automatic 
mode) an estimate of the full width at one- half of the maximum (FWHM) of a typical 
peak. 

The fitting routine does each peak in turn from the lowest to the highest energy. 
When two peaks are found to lie within four times the estimated FWHM, they are fitted 
simultaneously. Up to five peaks may be treated as a multiplet and done simultaneously. 
A more detailed description of the fitting procedure is given in appendix B. 

Following analysis of each line for energy and intensity, an identification is made 
by comparing the measured energy with a table of energies and atomic symbols. The 
elements which are included in this table are listed in table I. If a peak is identified as 
a line, a check is made to determine whether the Kg line can also be observed. If a 
line of the proper energy is located, the expected intensity of the line (determined 
from the intensity of the line) is subtracted from the candidate. The remainder 
of the peak is tested for other identification. No peak may be identified as an L series 
X-ray unless both the L a and lines are observed. 

The experimental intensity is converted to an areal density by using the stored effi- 
ciency function F(E). This function is assumed (for K series X-rays) to have the form 



log 


(6) 


7 

' i ~L V p(i) 

_F(E)J / j E (I-1) 

The parameters P(I) are determined by measuring the X-ray yield for a number of 
samples (approximately 10 to 15) of known areal density and performing a least- squares 
fit of the function given in equation (6) to that data. The parameters P(I) vary with any 
change in detector or source geometry and must be redetermined each time such a 
change occurs. A typical set of parameters P(I) is listed in table II. For heavier ele- 
ments, which are measured by means of D series X-rays, tabulated efficiencies are 
used for those elements which the program can identify. 

The estimates of errors in the computed areal density are based on the number of 
counts in each peak, given by 


N(E) = I(E)<r(E)Vir (7) 

where I and a (as defined in appendix B) are the height and width, respectively, of the 
Gaussian distribution derived from the data. The areal density, given by equation (5), 
then becomes 




I(E)a(E)VffF(E) 


N, 


Rh 


Therefore, assuming the error in Nj^ to be small gives 


51(E) 

+ 

6ct(E) 

+ 

6F(E) 

1(E) 


ct(E) 


F(E) 


( 8 ) 


(9) 


or 


61(E) 


5ct(E) 


1 dF(E)* E 

1(E) 

T 

ct(E) 


F(E) dE 


(10) 


The errors in each of the quantities 1(E), cr(E), and E are derived in the usual way 
from the residuals of the least-squares fitting process. 

Computed areal densities may be modified by an arbitrary multiplicative factor G, 
which may be entered by means of thumb wheel THl. The value of G may be between 
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0 and 9.9999, and its effect is included in all areal densities and errors listed as output. 
The principal purpose of this factor is to permit the user to convert to a convenient sys- 
tem of units if so desired. If G is set equal to 1, the results are in units of micro- 
grams per square inch. 

Output may consist of either one or two pages of teletype output as well as a record 
of the results on magnetic tape. These options are selected through the control panel 
toggle switches (p. 1, T14 up; p. 2, T12 up; and tape, T13 up). Samples of the two 
pages of teletype output are shown in figures 5 and 6. More detailed information on the 
results of the least-squares fitting process may be obtained by replying M T" (true) when 
the program inquires, through the teletype console, "Do you wish extra output? Out- 
puts " A sample of this augmented output is shown in figure 7. It includes the height, 
width, and location, as well as the standard deviation of each, for each line fitted by 
the program. 

A flow chart for the entire program is shown in figure 8. A guide to the loading of 
the program into the PD P-15 is given in appendix C. Appendix D contains the details of 
the overlay structure necessary in order to fit the program into the available core. 
Appendix E contains a listing of the FORTRAN sections of the program. 


CONCLUDING REMARKS 

The program described in this report permits trace element analysis of environ- 
mental samples at a rate which would make meaningful monitoring surveys feasible. 
Specifically, data analysis requires between 3 and 12 minutes, depending on the com- 
plexity of the spectrum, while data acquisition requires 5 to 30 minutes, depending on 
the primary radiation source and associated hardware employed. The program works 
reliably in an automatic mode; it changes samples, takes data, locates peaks, measures 
intensities, and identifies elements. As a result samples may be processed without an 
operator in attendance, which makes around-the-clock analysis possible. In addition, 
when operated under the supervision of the user, the program is capable of treating a 
wide variety of nonroutine samples, such as those containing elements not automatically 
recognized and those containing a number of heavy elements and thus having extremely 
complex spectra. For these, the program permits the extraction of meaningful data with 
a minimum amount of effort on the part of the user. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, December 2, 1976, 

506-25. 
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APPENDIX A 


AUTOMATIC PEAK LOCATING PROCEDURE 


If the operator so desires, initial estimates of peak locations are determined by 
the program, using a method described by Mariscotti (ref. 3). The method is based on 
the construction of smoothed second differences of the experimentally measured energy 
spectrum. The smoothed second difference is defined by 


i+m Z+m 

Sjd.m). V • • • (Nj +1 - 2N; + Nj.j) (Al) 

j=i-m h=Z-iq 

z 

or equivalently 


S.(z,m) = Cy(z, m)Nj (A2) 

allj 

where the C.j are weighting factors derived by Mariscotti, given by 

i+m 


Cj.(z, m)= ^ C £j(z-l>m) (A3) 



Z=i-m 




^0 if 

|j-i| — 2 


C i -(0,m) = < 

1 if 

Ii-i| =1 

(A 4) 


^2 if 

j =i 



The values of z and m used here are those found by Mariscotti to be most suit- 
able, namely, 


z = 5 


_ (o. 6r - l) 

2 
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The resulting function S then vanishes for a linear spectrum and is similar to the 
second derivative of a Gaussian function if Nj is a Gaussian function. The program 
assumes that a peak exists in channel i if three conditions are present: First, 


|S.| > 1.8 |6 Si | (A5) 

where 5S^ is the standard deviation in the function Sj, approximated by Mariscotti as 



where SNj is the standard deviation in the data point Nj. The program assumes 

6Nj = (N. + B.) 1/2 

where EL is the amount of background previously subtracted from channel i of the 
spectrum. Second, 


Sj < 0 


(A 6) 


Third, 


S i-1< 


I Sj | 


> 


S i+1 


(A7) 


The code does not reject any peak on the basis of its width as suggested by Mariscotti. 
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APPENDIX B 


i 


LEAST- SQUARES FITTING 


Once initial estimates of peak locations and widths have been obtained, a least- 
squares fitting is carried out by using the Gaussian function 


f(x) = I exp 


(x - x Q ) 2 


(Bl) 

* 


to represent that portion of the spectrum being treated. In the event that several peaks 
fall within 4FWHM of each other, they are fitted simultaneously as the sum of as many 
as five Gaussian distributions, each of which has independently variable width. 

For initial estimates, Xq is taken to be the channel i either found by the automatic 
search routine described in appendix B or supplied through the light pen. The initial 
estimate of the intensity I is taken to be the number of counts in that i** 1 channel, and 
a is calculated from the FWHM input (for a Gaussian distribution, a = FWHM/1. 665). 

Fitting is carried out by standard methods. In the event a satisfactory fit cannot 
be obtained within a specified number of iterations, the peak of highest energy is 
dropped from the multiplet, and the fitting process is restarted for the remainder of 
the peaks. In the event that a single peak cannot be satisfactorily fitted, a message to 
that effect is included on output page 1, and the peak is assigned intensity zero and a 
location equal to the starting channel number. 


d'i 
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APPENDIX C 


LOADING THE PROGRAM 

As presented in this report the program (entitled SAVLOC) is too large to be loaded 
directly into a 16K PDP-15 computer. As a result it can only be used as a series of 
overlays, constructed by using the routine CHAIN (ref. 4). A listing of the execute file, 
titled SAVLOC, is given in appendix D; the listing includes the resident code as well as 
the structure of each of the links and a machine constructed map of the file . 

At execution time the user must assign handlers to three slots in the device assign- 
ment table (dat): dat slot 1, for writing of raw data; dat slot 2, for reading of raw data; 
and dat slot 3, for writing of results of analysis (fig. 6). Dat slots 1 and 2 are normally 
assigned the same magnetic tape (dectape) unit, and dat slot 3 is assigned a second 
dectape unit. 

An example of the entire loading sequence is as follows: 

KMS15 V5C 

_$A DTE 6 1, 2/DTE7 3 

_$E SAVLOC 

EXECUTE V4A 
FNAME(A) = 


where the underlined expressions are output from the computer. 
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APPENDIX D 


OVERLAY STRUCTURE AND MEMORY MAP FOR EXECUTE FILE SAVLOC 


\CHAIN 


CHAIN V5B 

NAME XCT FILE 
>SA VLOC 

LIST OPTIONS 8 , PARAMETERS 

>PGR, 16K, SZ 

DEFINE RESIDENT CODE 

>SI MUL2 , CAMADC, IPK3, PICKER 

DESCRIBE LINKS & STRUCTURE 

>LK1 =PEDE, RITE 

>LK2 = FIT ALL, FITAL1 , GAUSS 

>LK3 = ECAL1 , ECAL2 

>LK4=PEG0UT, XKYLD 

>LK5=LAREL, SOFT 

>LK6 = S, SCAL1 , LOCUM 

>LK1 : LK2: L K3 : LK4 : LK5: LK6: PEED: MAKER: SET1 : StJBTR: OUTER: 
- MULT: SIT2: EPPO: ENERG: I PENT 


LINK TABLE 

37225-37636 004 P 

RESIDENT CODE 
SI MULE 36527-37224 00476 
CAMADC 35731-36526 00576 
IPX3 34740-35730 00771 
PICKER 34305-34737 00433 
RELAON 34226-34304 00057 
RTIMER 34146-34225 00760 
RFCV16 34 1 26-34 1 45 00020 
SEND16 34 1 l 3-34 1 25 00013 
CAMAC. 33631 -341 1 2 00262 
PANEL 33477-33630 0013° 
CROSS 33347-33476 00130 
HI ST 01 33226-33346 00121 
PE NSW 33200-33225 00026 
VPVECT 33014-33177 00164 
VIDEO. 32515-33013 00277 
RT I ME 32470-32514 00025 
CLOCK. 32372-32467 00076 
R I TEST 32350-32371 00022 
3ETRLK 32267-32347 00061 
IROOL 32224-32266 00043 
ZERO 32203-32223 00024 
.SS 32130-32177 00050 
FLOAT 32117-32127 30011 
FL0AT2 32071-32116 00726 
.BR 32024-32070 03045 
FIOPS 31302-32023 00522 
INTEAE 31 171-31301 001 1 1 
PELEAE 30057-31170 01112 
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OTSEP 277 1 2-27777 00066 
.DA 3P5H21 -30056 00036 
D I PPL 300 1 6-30022 00003 
EPP 30012-30015 00004 
STUFF 27464-2771 1 02226 
nlJP 27445-27463 00017 
FINOUT 27345-27444 20100 

LINK — LK1 

REDE 27277-27344 00046 
RITE 27135-27276 00142 
RELAOF 27061-27134 00054 
B CD EXT 27012-27060 00047 
INPUT1 2672.3-27011 00067 
XPEAD 26407-26702 00314 
PACKER 24306-26406 00101 
MOD 26265-26305 00021 
BCDIO 23630-26264 02435 
BIN 10 23403-23627 00225 
FILE 23124-23402 00257 
.CB 23105-23123 00017 

LINK — LK2 

FITALL 26136-27344 01207 
FITAL1 26005-26135 00131 
GAUSS 25277-26004 00506 
MATINV 24736-25276 00341 
GENRLK 24663-24735 00053 
LOC 24654-24662 00007 
ABS 24635-24653 00017 
RMNMX 24511-24634 00124 
.RC 24415-24510 00074 
EXP 24377-24414 00016 
. EF 24272-24376 00105 
.EC 24220-24271 00052 
.R DA 24130-24217 00070 
FITALC 24127-24127 00001 

LINK — LK3 

ECAL1 27001-27344 00344 
ECAL2 26446-27000 00333 
INPUT1 26357-26445 00067 
XREAD 26043-26356 00314 
PACKER 25742-26042 00101 
BCDIO 23305-25741 02435 
.CB 23266-23304 00017 

LINK — LK4 

REGOUT 26470-27344 00655 
XKYLD 26343-26467 00125 
.BC 26247-26342 00074 
.BE 26237-26246 00010 
.PR 26201-26236 00036 
SQRT 26066-26200 00113 
.EE 25703-26065 00163 
.EF 25576-25702 00105 
.EC 25524-25575 00052 
.RDA 25434-25523 00070 


IL 


BCDIO 

.CB 


22777-25433 02435 
22760-22776 00017 


LINK — LK5 

LABEL 24757-27344 02366 
SORT 24621-24756 00136 
ABS 24602-24620 00017 


l 


LINK 
S 

SCAL1 
LOCUM 
IFIX 
ABS 
LIMIT 
.BC 
SQPT 
DIFF 

LINK -- REED 

PEED 26775-27344 00350 
BCDEXT 26726-26774 00047 
INPUT1 26637-26725 00067 
XPEAD 26323-26636 00314 
PACKER 24222-26322 00101 
MOD 26201-26221 00021 
BCDIO 23544-26200 02435 
BINIO 23317-23543 00225 
FILE 23040-23316 00257 
.CB 23021-23037 00017 

LINK — MAXER 

MAKER 27275-27344 00050 


— LK6 

27117-27344 00226 
24651-27116 00246 
24730-26650 01721 
24712-24727 00016 
24673-24711 00017 
24650-24672 00023 
24554-24647 00074 
24441-24553 00113 
24220-24440 00221 


LINK — SET 1 

SFT 1 27003-27344 00342 


LINK -- SUBTP 
SUBTP 26527-27344 00616 
INPUT1 26440-26526 00067 
XPEAD 26124-26437 00314 
PACKER 26023-26123 00101 
BCDIO 23366-26022 02435 
.CR 23347-23365 00017 


LINK -- OUTER 
OUTER 27041-27344 00304 
BCDIO 24404-27040 02435 
.CB 24365-24403 00017 

LINK — MULT 

MULT 26665-27344 00460 
IFIX 26647-26664 00016 

LINK — SIT2 

SIT2 26032-27344 01313 
BCDIO 23375-26031 02435 
.CP 23356-23374 00017 


i 
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LINK - 

- ERRO 


ERRO 

27033-27344 

00312 

BCDIO 

24376-27032 

02435 

.CB 

24357-24375 

0001 7 

LINK - 

- ENERG 


ENERG 

26500-27344 

00645 

LINK - 

- IDENT 


I DENT 

26066-27344 

01 25 7 

RCDEXT 

2601 7-26065 

00047 

PACKER 

25716-2601 4 

00101 

APR 

25677-25715 

0001 7 

MOD 

25656-25676 

00021 

BCDIO 

23221-25655 

02435 

file 

22742-23220 

00257 

.CB 

22723-22741 

00017 

BLANK 

COMMON 


.XX 

12723-22722 

1 0000 

CORE REQ' D 



12723-37636 

24714 

KMS1 5 

V5C 
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APPENDIX E 


SOURCE LISTING 


C PROGRAM TO SIMULTANEOUSLY TAKE AND ANALYZE DATA 

INTEGER DATAD( 4096) , GATEXC64 ,7) , O.RGN(8) , DATUM( 1 024) 

REAL FNAME(2) , LT, LT1 
LOGICAL ON, EXIST, BUTTNS,B ITEST 
EXTERNAL ADCON 
COMMON DATAD 

EQUIVALENCE ( GATEX( 1 , 1 ) ,DATAD( 3573) ) , (DATAD( 1 ) , DATUM ( I ) ) 
DATA ND,NLAST/1 024,1/, 

1 FNA ME C 1 ) , FN AME ( 2) /5 HSPECT , 1 H / 

CALL CLKON 
CALL PSTART 

2 CALL REDE( FNAME, NEXT, NLA ST) 

3 I4=IROTOR(4) 

I3=IR0T0R(3) 

I2=IR0T0R(2) 

EXIST=.TRUE. 

NY=IR0T0R(4)+6 
NF=256*I2-255 
NX: 13+7 

CALL ZERO(DATUM) 

CALL CLKON 
CALL VP STOP 

CALL HISTOl (DATUM, NF, NX, NY, NSAW, 1 ) 

CALL VPSTRT 
CALL ECAMCK 
CALL RTCOFF 
IF(ON) CALL ADCOFF 
104 MODE=0 

IRT= ITHUMBC2) 

ILT = 0 

IF( NEXT.GT. NLAST) GO TO 100 

CALL ADCSTRC DATAD, GATEX, ORGN, ILT, ON, MODE) 

CALL RTCSTPCIRT, ADCOFF) 

CALL RTIME(TS) 

CALL LAMPON(l) 

100 CALL BUTTNS(IB) 

IFCI3.EQ.1) GO TO 101 

IF( IR.EQ.2) GO TO 501 

IF( . NOT. ON) GO TO 96 

102 IF(IB.EQ,4) GO TO 98 

CALL TOGGLS(IT) 

IF(BITEST(IT, 15) .AND. EXIST) GO TO 101 
GO TO 100 

101 CALL VP STOP 
CALL IPK2CEXIST) 

CALL VP STOP 

CALL HISTOl (DATUM, NF, NX, NY, NSAV, 1 ) 

CALL VPSTRT 

99 I F ( ON) 30 TO 100 

96 IF( NEXT.GT. NLAST) GO TO 102 

98 CALL RTCOFF 

CALL ADCOFF 
CALL LAMPOF(l) 
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CALL TOGGLSC IT) 

IFCIROTORCl) .E0.2) CALL RELAON (3) 
CALL RI TEC FNAME, NEXT, TS, IRT, ILT, ND) 
NEXT: NEXT+1 
CALL TOGGLSC I T) 

IFCBITESTCIT,15)) GO TO 3 
CALL LAMP0NC2) 

500 CALL BUTTNSCIB) 

IFCIB.NE.2) GO TO 500 

502 CALL LAMP0FC2) 

501 CALL T03QLS< IT) 

IF(BITESTCIT,17)) GO TO 2 
GO TO 3 

END 


SUBROUTINE IPK2C EXIST) 

C SUBROUTINE TO DO PEAK ANALYSIS AFTER LIGHT PEN SELECTION, 

C FOR USE WITH SIMUL2, TO TAKE DATA SIMULTANEOUSLY, 

C USED TO BIULD XCT FILE 

INTEGER LI STC520) , DATA (1024) ,BKGND(1024) , GATE( 5 , 30 ) , CH0 (30 ) , 

1 CHMAX, DATADC4096) ,BEXT, IRC5) 

REAL FNAME (2) , PARAMO 5) , DP( 1 5) , POI NT( 1 00) ,B NAME(2 ) ,XCH(2) , 

1 ELEC 30) ,ENORM(30) ,BNORMX(30) ,DR( 30) 

LOGICAL EXIST, DOUBLE, UNDOUB, TRIPLE, UNTRIP, QUAD, UNQUAD, QUIN, 

1 UNQUIN, BUTTMS, PENSW, OUTPT, BITEST 

EXTERNAL GAUSS 
COMMON/D I SPL/NF, NX, NY 
C0MM0N/ERR/PHI2, IMI N, I TELL 
COMMON/STUFF/GATE 

COMMON/MUP/DOUBLE, UNDOUB, TRIPLE, UNTR IP, QUAD, UNQUAD, RUI N , UNOUI N , 
1 JPEA K, GAM, IDUB, ITRIP, IQUAD, IQUIN 

COMMON/FINOUT/ELE, A1 , A2 
COMMON DATAD 

EQUIVALENCE (DATA (I ), DATA D( 1 025 )),( RKGN D( 1 ) ,DATAD(2049) ) , 

1 (LISTC1 ) , PAT ADC 3073) ) , (CH0 ( 1 ) , DATADC 4021 ) ) , 

2 (POINTd ) , DATADC 3593) ) , CENORMC1 ) , DATA DC 3800) ) , 

3 CBNOPMXC1 ) , DATADC 3861 )), C DB ( 1 ), DATADC 392 1 ) ) , 

4 C IRC 1 ), DATADC 4051 )), (PARA MCI), DATADC 4056) ) 

CALL LAMPONC3) 

CALL REEDCl , FNAME, NEXT, EXIST) 

IFC EXIST) GO TO 3 
CALL LAMPOFC3) 

RETURN 

3 NEXT: NEXT+1 

2 CALL MAXERC DATA, CH MAX) 

CALL SET1 (CHMAX, DATA, POINT) 

ITELL=1 

CALL FI TALLC GAUSS, POINT, PAR AM, DP, Q, LIST) 

CALL OUTERC PARAM, LI ST, A NORM, A p NORM, CHMA X) 

CALL REEDC2,B NAME, BEXT, EXIST) 

CALL MAXEPCBKGNP, CHMAX) 

CALL SET1 (CHMAX, BKGND, POINT) 

CALL FITALLCGAUSS,POINT,PARAM,DP,Q, LI ST) 

I TELL: 0 

CALL SUBTRC A NORM, A R NORM, CHMAX, OUTPT, GAM) 

CALL PICKER 
CALL TOGGLSC IT) 

IFC PI TEST (IT, 10)) CALL VPSTOP 
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IFCBI TEST< I T, 1 1 )) CALL LOCUM 
CALL LAMPOF(7) 

IFCBITESTCIT, 16)) GO TO R7 
DO 88 J=1 ,2 
CALL ECALl(J) 

CALL FITALLC GA'JSR, POI NT, PARAM, DP, Q, LI ST) 

XCH( J)=FLOAT(CH0( J) )+PARAM(2) -4. 

88 CONTINUE 

87 CALL ECA L2 ( XCH , A 1 , 42 ) 

C BEGIN FITTING PEAKS 

80 CALL LA MPOF (7) 

DO R9 J=3, JPEAK 

81 CALL MULT(IMAX,IMIN,J) 

CALL SIT2( IMAX,IMIN, J.ANORM) 

CALL FI TALL( GAUSS, POI NT, PARAM, DP, Q, LI ST) 
NPARAM=LIST(6) 

DO R5 1=1 , NPARAM,3 

IFCPARAMCI+U .GT. FLOAT ( IMAX-IMIN+5) ) LISTC4) =7 
IF(PARAM(I+1 ).LT.-5.) LIST(4)=R 

85 CONTINUE 

IF(LIST(4) .EO.0) GO TO 86 
CALL ERROC J, IMIN,DP) 

IFC. NOT. DOUBLE) GO TO 86 
GO TO 81 

86 CALL ENERG( I MI N, 41 ,A2 , J, PP) 

CALL LABEL(J, JPEAK) 

CALL REGOUTC OUTPT, ANORM, IMIN,DP,Q, J) 

89 CONTINUE 

CALL IDENTC ELE , ENORM,B NORMX, DB, JPEAK, FNAME, NEXT) 
CALL LAMPOFC3) 

CALL VPSTRT 

RETURN 

END 


subroutine PICKER 

INTEGER GATE (5, 30) ,CH0 (30) ,DATA( 1 024) , DATADC4096) , OFF, X, Y, 

1 CHMAX 

LOGICAL DOUBLE, UNDOUB, TRIPLE, UNTR IP, QUAD, UNQUAD, QUIN, UNQUIN, 

1 BUTTNS,PENSW, BITEST 

COMMON/STUFF/GATE 

COMMON/MUP/DOUBLE, UNDOUB, TRIPLE, UNTRIP, QUAD, UNQUAD, QUIN, 

1 UNQUIN, JPEAK, GAM, IDUB, I TRIP, I QUAD, I QUIN 

COMMON/DI SPL/NF, NX, NY 
COMMON DATAD 

EQUIVALENCE (DATAC1 ) ,DATAD(1 025) ) , (CH0C1) ,DATAD(402i ) ) 

DATA OFF/-1 / 

100 FORMATd X,2H0K) 

ICHK=0 
6 IS1 = NX 

I S2 = NY 
IS3=NF 

CALL SCAL1 ( NX, NY, NF) 

CALL TOGGLS(IT) 

IF( I CHK.NE, 1 ) GO TO 66 

IFCCISl .NE.NX) .OR.CIS2.NE.NY) .OR. (IS3. NE. NF) ) GO TO 66 
GO TO 65 

66 CALL VP ST OP 

CALL HISTOl ( DATA , NF, NX, NY, NHIT, 1 ) 


18 



CALL VPVECTCGATE,30,IHIT,2) 
CALL VPSTRT 

65 IF(8ITESTCIT,16)) RETURN 

IF(ICHK.EQ.l) GO TO 12 
IFC ICHK.EQ.-l ) GO TO 13 

7 CALL ZEROCGATE) 

CALL ZERO(CH0) 

JPEAK=0 

CALL CR0SS(X,Y,0FF,3) 

8 JPEAK= JPEAK+1 

12 IFCPENSW(.TRUE.)> GO TO 9 

ICHK=1 

GO TO 6 

9 GATE Cl , JPEAK) = NHIT 
CH0C JPEAK) =NH1T 
GATEC2, JPEAK) =DATAC NHIT) 
GATEC3, JPEAK) =0 
GATEC4, JPEAK) =128 

GATE C5, JPEAK) =-10 

ICHK=-1 

GO TO 6 

13 CALL LAMPONC 7) 

10 CALL BUTTNSC IB) 

IFCIB.EQ.128) GO TO 1 1 
IFCIB.EQ.32) GO TO 14 
IFCI3.E0.64) RETURN 

GO TO 10 

14 JPEAK= JPEAK-1 

1 1 CALL LAMPOFC7) 

GO TO 8 
END 


SUBPOUTI NE REDEC FNAME, NEXT, NLA ST) 
REAL FNAME (2) 

CALL I NPUT1 ( FNAMEC 1 ) ,5HFNAME,6) 
CALL I NPUT1 C NEXT, 4HNEXT, 1 ) 

CALL I NPUT1 ( NLAST,5HNLAST , 1 ) 

RETURN 

END 


FUNCTION SCICHNL) 

INTEGER DA TA ( 1 024) , PAT AD (40 96) 

REAL CC36) ,BC36) 

COMMON/DIFF/C,B, M 
COMMON DA TAD 

EQUIVALENCE C DATAC 1 ) , DATA DC 1 025) ) 

Ml =5* W-2 

X=C( 1 )* FLOAT C DATAC ICHNL) ) 

DO 41 J = 1 , Ml 
II = ICHNL+J 
12 = I CH NL- J 
J1 =J+1 

IFCI1 ,GT, 1024) GO TO 42 
IFCI2.LT. 0) GO TO 42 

XrX+CC J1 ) * FLOATC DATAC 1 1 ) )+CC J1 )* FLO A T( DATAC 12) ) 
GO TO 41 
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42 X=0 

41 CONTINUE 

S= X 

RETURN 


SUBROUTI NE RITE (FNAME, NEXT, TS, IRT, ILT, ND) 

REAL FNAMEC2) 

INTEGER DATADC4096) 

COMMON DATAD 

CALL BCDEXTC FNAME, NEXT) 

CALL ENTER(2, FNAME) t 

WRITE (2) TS, IRT, ILT, ND, < DATADC I) , 1 = 1 , ND) 

CALL RELA0FC3) 

WRITEC6, 106) FNAME, TS, IRT, ILT 
106 FORMATC1 X,2A5,5X, FI 0.5, 16 ,3X, 16) 

CALL CL0SE(2) 

RETURN 

END 


SUBROUTINE GAUSSCR, X, P,DP, V,KK, J) 

SUBROUTINE FOR US EVITH FITALL TO FIT UP TO 4 
GAUSSIAN PAKS WITH OR WITHOUT BACKGROUND 
REAL X(10e),P(15),DP(15) 

INTEGER BKGNDC1024) ,DATA( 1 024) , DATAD( 4096) 
C0MM0N/ERR/PHI2 , 1 LOWER, I TELL 
COMMON DATAD 

EQUIVALENCE ( DATAC 1 ), DATADC 1 025) ), (RKGNDC 1 ), DATADC2049) ) 
NP=KK/3 

IFCKK-3* NP.EQ.0) GO TO 1 
R = XC J)-PC1 )* FLOATC J)-P(2) 

DP( 1 ) =-FLOA TC J) 

DP(2)=-1 . 

N = 3 

GO TO 3 

1 N=1 

R = XC J) 

3 DO 2 I=N,KK,3 

TP=P(I)* EXP (-( FLOATC J) -PC 1+1 ) )**2/P( I+2)**2) 

R=P-TP 

DP( I ) =-TP/P ( I ) 

DPCI+1 ) =- TP*2 .* ( FLOAT ( J)-P(I+l))/P(I+2)**2 
DP(I+2)=DP(I+1)*( FLOATC J)-P(I+1 ))/P(I+2) 

2 CONTINUE 

IERR= ILOW ER+J-3 

I F ( I TELL. EO .0) V= C FLOATC BKGND( IERR)+1 ) ) 

IF< ITELL.EO.l ) V=l. 

RETURN 

END 


SUBROUTINE ECALl(J) 

C SUBROUTINE CALL ED BY IPK2 TO DETERMINE ENERGY CALIBRATION. 

INTEGER CH0C30) , LI ST( 520) , DATADC 4096) , DA TA< 1024) 

LOGICAL BITEST 

REAL PARAMO 5) ,DPC 1 5) ,POINT( 1 00) , XC2) 

COMMON DATAD 

COMMON /ERR/PH 12, 1 MIN, I TELL 
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EQUIVALENCE ( DATA C 1 ) , DATAD( 1 025) ) , C LI STC ! ) , DATA r >C30 73) ) , 

1 ( CHQ ( 1 ) ,DATADC4021) ) ,CPOI NTH ) ,DATADC 3593)) , CPARAMC1) 

1 .DATADC4056)) 

CALL TOGGLSCIT) 

IFCBITESTCIT, IS)) RETURN 
LISTCl )=7 
LI STC2) =3 
LI STC 3) = 1 5 
LISTC4) -?■ 

LIST<5)=3 
LISTC6) =3 
IK=CH0C J) 

PARAMCl ) =FLOA TC DATAC I K) ) 

PARAMC2) =4, 

PA PA II (3) =4.5/1 .S65 

I MI N = CH0 ( J) -3 

DO 3 1=1,7 

IC=CHP( J) -4+1 

POINTC I) = FLOAT ( DA TA( IC) ) 

3 CONTINUE 

RETURN 
END 


SUBROUTINE ECAL2CX,A1 , A2) 

C SUBROUTINE TO DO ENERGY CALIBRAION WITH IPK2 

REAL XC 2) 

LOGICAL BITEST 
CALL TOGGLS(IT) 

IFCBITESTCIT,1S)) GO TO 99 
WRITE (6,100) 

CALL INPUT1 < El ,2HE1 ,2) 

CALL INPUT1 CE2,2HE2,2) 

A1 = CE1-E2)/CXC1)-XC2)) 

A2=E1 -A 1 *XC 1 ) 

99 G=FLOAT( I THUMB Cl ) ) /I 0000 . 

WRI TECS , 1 03) G 

103 FORMA TCI X,2HG= , FS.4) 

IFCRITESTCIT, 14)) WRITECS,101) 

IFCRITESTCIT, 14)) WRITECS,102) 

RETURN 

101 FOPMATC ///, SH PEAK, 9H CHANNEL, 12H STRENGTH, 

1 10H ENERGY, 1SH CONCENTRATION,! 0H ERROR, 

p pj-j flem ) 

102 F0RMATC39X,f4HC UGRAMS/ 1 N**2) ) 

100 FORMATC 1 X,25HENERGIES FOR CALIBRATION:) 

END 


SUBROUTINE REGOUTCOUTPT, ANORM, I MIN, DP, (3, J) 

* INTEGER LISTC520) .DATADC409S) 

REAL PARAMO 5) , ENORMC30) , ELEC30) ,DP< 1 5) .BNORMXC30) ,DBC30) 

LOGICAL OUTPT, BITEST 

COMMON DATAD 

COMMON/FI NOUT/ELE,Al ,A2 

EQUIVALENCE C LISTCl ) ,DATADC3073) ) , CENORMCl ), DATA DC 3 800) ) , 

1 CBNORMXC1 ) ,DATADC3B61 )) , CDRC 1 ) , DATADC3921 )) , 

2 CPARAMC1 ) ,DATADC405S)) 

NPARAM=LI STC6) 
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DO 92 1=1 ,NPARAM 
DP( I) =SQRT( (Q)*DP(I)> 

92 CONTINUE 

DO 91 I=1,NPARAM,3 
IP2=I+1 

AN0RMX=PARAMCI)*PARAMCI+2)*1 ,7725/ANORM 

BNORMXC J) = ANORMX*XKYLDC ENORMC J) )* FLOAT C ITHUMBC 1 ) )/l 0000 . 
DBCJ)=CDPCI)/PARAMCI))+CDPCI+2)/PARAMCI+2) ) 
DY=A1*DPCIP2)/XKYLDCEN0RMCJ)> 

DY=DY*CXKYLDCENORMCJ))-XKYLDCENORMCJ)+0.02))*50. 

DB C J) = C DB C J) +DY) *B NORMXC J) 

CALL TOGGLSC IT) 

IFCBITESTCIT,14)> WRITEC6,125) J, PARAMC I P2) ,ANORMX, ENORMC J) , 
1 BNORMXC J) ,DBCJ) ,ELECJ) 

J= J+l 


91 CONTINUE 

J=J-1 

IFCOUTPT) WRITEC6, 1 09) Cl, PARAMC I) , DP Cl), 1=1, NPARAM) 

109 F0RMATC1X,2HPC,I2,2H)=,F12.3,5X,6HDELTA=,F12.3) 

1 25 FORMATC1 X, 14 ,2 X, FB. 3, FI 1 . 6 ,3X, FB. 3,3X, El 1 . 5,2X,E1 1 . 5,3 X, 1 A5) 

RETURN 
END 

C FUNCTION XKYLD, TO TRANSFORM RATIO TO RHODIUM TO 


1 NANO GRAMS/ 0(^*2 

FUNCTION XKYLDCE) 

REAL PCS) 

DATA PC1),PC2),PC3),PC4),PC5),PC6),PC7)/ 

1 1 . 2347, 30. R51 6, -306.138,1 980.348,-6477.202, 

2 10661 .801 ,-7106. *53/ 

1 AMT=0. 

DO 2 1=1,7 

AMT=AMT+PC I)/CE)**CI-1 ) 

2 CONTINUE 
XKYLD=10.**AMT 
RETURN 

END 


SUBROUTINE LA8ELC J, JPEA K) 

REAL ENORMC30) , ELEC30) , X, MN, CRMN, N I, N I BET, MO, BNORMXC 30) , D3C30) 
INTEGER CH0C30) , LI STC 520 ) , DATADC 4096) 

LOGICAL CHECK 
COMMON/FINOUT/ELE, A1 ,A2 
COMMON DA TAD 

EQUIVALENCE CCH0C1 ) , DATADC 402 1 )), C LI STC 1 ), DATADC 3073) ) , 

1 C ENORMC 1 ) , DATADC 3800) ) , CBNORMXC 1 ) ,DATADC386 1 ) ) , 

2 CDBC 1 ) , DATADC 3921 ) ) 

DATA CL,AR,K,SC,CA,TI/2HCL,2HAR,1 HK,2HSC,2HCA,2HTI/, 

1 V, CR, MN, FE, CO, NI/1 HV,2HCR,2HMN,0HFE,2HCO,2HN I/, 

2 CU,Z N, TLL, PBL, SE,BR/2HCU,2HZN,3HTLL, 3HPRL,2HSE,2HBP/ , 

3 THL,RB, SR, ZR, MO/3 HTML, 2 HRB, 2HSR,2HZR, 2HM0/ , 

4 CASC, SCTI , TI V, VCR, CRMN/5 HCA- SC ,5 HSC- TI , 4 HTI-V, 

5 4HV-CF,5HCR- MN/, 

6 FEBET, NIBET, CUBET/5HFERET, 5 HN I BET , 5HCUPET/ , 

7 ZNBET,BRBET, UNKNW/5HZNPET , 5 HBRBET , 5 H / 

NELE=CLISTC6)/3) + J 

DO 2 I = J, NELE 
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E= ENORMCI) 

ELEC I) =UNKNW 

IF ( ABSCE-2 .62) .LE.0.06) ELE(I)=CL 
IFC A3SCE-2. 95). LE.0.06) ELECI) = AR 
IFCABSCE-3. 31) .LE.0.06) ELECI)=K 
20 IFCABSCE-3. 69) .LE.0.06) ELECI)=CA 

IFCABSCE-4. 09) .LE.0.06) CALL SORTCCHECK, 3. 69, JPEAK) 
IFC.NOT. CHECK. AND. CABSCE-4. 09). LT. 0.06)) ELECI)=SC 
IFC CHECK) ELECI) =CASC 
CHECK= .FALSE. 

IFC ABSC E- 4^51). LE. 0.06) CALL SORTCCHECK, 4.09, JPEAK) 
IFC. NOT. CHECK. AND. CABSC E-4.51 ) .LT.0. 06)) ELECI) =TI 
IFC CHECK) ELEC I) = SCTI 
CHECK=. FALSE. 

IFCABSCE-4. 95) .LE.0.06) CALL SORTCCHECK, 4.51 , JPEAK) 
IFC. NOT. CHECK. AND. CABSCE-4. 95) .LE.0.06) ) ELECI) = V 
IF C CHECK) ELECI) =TIV 
CHECK=. FALSE. 

IFCABSCE-5. 41). LE.0.06) CALL SORTCCHECK, 4. 95, JPEAK) 
IFC. NOT. CHECK. AND. CABSC E-5.41) .LE.0.06)) ELECI) =CR 
IFC CHECK) ELECI) =VCR 
CHECK=. FALSE. 

IFCABSCE-5. 900) .LE.0.06) CALL SORTC CHECK, 5 . 4 1 , JPEA K) 
IFC .NOT. CHECK. AND. CABSC E- 5. 90) .LE.0.06) ) ELECI) =MN 
IFC CHECK) ELECI):CRMN 
CHECK=. FALSE. 

28 IFCABSCE-6. 40) .LE.0.06) ELE Cl) =FE 

IFC ABSCE-6 . 93) .LE.0.06) ELECI) =CO 
IFCABSCE- 7. 06) .LE.0.06) ELECI) =FERET 
IFCABSCE-7.48) .LE.0.06) ELECI)=NI 
IFCABSCE-8. 05) .LE.0.06) ELECI)=CU 

IFCABSCE-S. 26) .LE.0.06) .CALL SORTC CHECK, 7 .48 , JPEAK) 
IFC CHECK) ELECI) = NIBET 
CHECK:. FALSE. 

IFCABSCE-8. 63) .LE.0.06) ELECI) =7N 

IFCABSCE-R. 90). LE.0.06) CALL SORTCCHECK, 8. 05, JPEAK) 
IFC CHECK) ELECI) =CUBET 
CHECK:. FALSE. 

IFCABSCE-9. 57) .LE.0.06) CALL SORTCCHECK, 8. 63, JPEAK) 
IFC CHECK) ELE Cl) N3ET 
CHECK: .FALSE. 

IFCABSCE-10. 25) .LE.0.06) CALL SORTC CHECK, 1 2 .2 1 , JPEA K) 
IFC CHECK) ELECI) =TLL 
CHECK=. FALSE. 

IFCABSCE-10. 55) .LE.0.06) CALL SORTC CHEC K, 1 2 . 6 1 , JPEA K) 
IFC CHECK) ELECI) :PBL 
CHECK:. FALSE. 

34 IFCASSCE-11 .22) .LE.0.06) ELECI)=SE 

IFC ABSCE-1 1 .89) .LE.0.06) ELEC I) =BR 

IFCABSCE-I2. 21 ) .LE.0.06) CALL SORTC CHECK, 1 0 .27 , JPEA K) 
IFC CHECK) ELECI) :TLL 
CHECK:. FALSE. 

IFCABSCE-12. 61) .LE.0.06) CALL SORTC CHECK, 1 0 . 55 , JPEA K) 
IFC CHECK) ELECI) :PBL 
CHECK:. FALSE. 

IFCABSCE-12. 90) .LE. 0.10) CALL SORTC CHECK, 16 .20, JPEA K) 
IFC CHECK) ELECI) =THL 
CHECK:. FALSE. 

IFCABSCE-1 3.29) .LE.0.06) CALL SORTC CHECK, 1 1 . 92, JPEA K) 
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I F < CHECK) elfci)=rrbet 
CHECK: .FALSE. 

IFCARSCE-13.39) .LE.0.06) ELECI)=RB 
IFC ARSCE- 1 4 . 1 6) .LE.0.06) ELECI)=SR 
IFCARSCE- 15.77) .LE.0.06) ELECI)=ZR 

IFC ARSC E-l 6 .20) . LE.0 • 1 0) CALL SORTC CHECK, 12.90, JPEA K) 
IF( CHECK) ELE(I) =THL 
CHECK: .FALSE. 

42 IFCARSCE-17.47) .LE.0.06) ELECI)=MO 

2 CONTINUE 

RETURN 
END 


SUBROUTINE SORTCCHECK, X,.I?E4K) 

INTEGER CH0 (3?) , DATADC409S) 

LOGICAL CHECK 

REAL ELEC30>,RNORMXC30>,DBC30) 

COMMON/FI NOUT/ELE, A1 , A2 
COMMON DA TAD 

EQUIVALENCE CCH0 Cl ) , DATADC4021 ) ) , (BNORMXC I ) , DATADC3861 ) ) , 
1 CDB ( 1 ), DATA DC 3921 ) ) 

CHECK:. FALSE. 

DO 2 1=1, JPEA K 
E: A 1 * FLOA TC CH0C I) ) +A2 
IFCABSCE-X) .LE.0.10) GO TO 5 
GO TO 2 
5 I: JPEA K 

CHECK:, TRUE. 

2 CONTINUE 

RETURN 
END 


SUBROUTI NE SCAL1 C NX, NY, NF) 

SUBROUTINE TO CHANGE SCALES ON DISPLAY, 

FOR USE WITH IPICK 

INTEGER GATE C5 ,30) ,DATAC 1 024 ) , CH0 C30 ) , DATADC 4096) , CHMAX 

COMMON/STUFF/GATE 

COMMON DATAD 

EQUIVALENCE C DATAC 1 ) , DATAD C 1 025) ) , ( CH0 C1),DATADC 402 1 ) ) 
I2=IROTORC2) 

NF=LIMITC1 ,256*12-255,1024) 

I3=IR0T0PC3) 

NX=LIMITCB, 13+7, 1 2) 

I4=IR0T0RC4) 

NY=LIMITCP, 14+6, IP) 

DO 10 1=1 ,30 
LI =CH0C I) 

IFCL1.EQ.0) LI = 1 

GATE Cl , I) =LI MITC0.CCH0CI )-NF)*B/2**I3, 1023) 

GATEC2, I)=DATACL1 )/C2**I4)*l6 
10 CONTINUE 

RETURN 
END 
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SUBROUTINE REEDCI TYPE, FNAME, NEXT, EXIST) 

.SUBROUTINE TO READ DATA FROM TAPE FOR USE WITH IPK2 
AND SIMUL2 

INTEGER DATAC1024) , DATADC 4096), BKGNDC 1024) 

REAL FNAMEC2) 

LOGICAL EXIST,B I TEST 
COMMON DATAD 

EQUIVALENCE C DATA C 1 ), DATAD ( 1 025) ) , CB KGNDC 1 ) , DATADC 2049) ) 
CALL TOG GL SC IT) 

IFCBITESTCIT, 15)) GO TO 3 

1 CALL INPUT1 CFNAMEC1 ) ,5HDNAME,6) 

CALL INPUT1 C NEXT,4HNEXT, 1 ) 

3 CALL BCDEXTC FNAME, NEXT) 

CALL FSTATC1 , FNAME, EXIST) 

I FC EXIST) GO TO 2 
WRITEC6, 100) FNAME 
IFCBITESTC IT, 1 5) ) RETURN 
GO TO 1 

2 IFCITYPE.EO.l) CALL ZEROCDATA) 

IFCITYPE.EO.l) WRITEC6 , 1 01 ) FNAME 

101 FORMATC 1 HI , 1 2H PUN NUMBER ,2A5) 

IFCITYPE.EQ.2) CALL ZEROCBKGND) 

CALL SEEKC1 , FNAME) 

IFCITYPE.EO.l )READC1 ) TS,IRT,ILT, ND, CDATAC I) ,1=1,1024) 
IFCITYPE.EQ.2) READC1) TS, IRT, ILT, ND, CBXGNDC I) , I =1 , l 024) 
RETURN 

100 FORMATC 1 X,6H0 FI LE ,2A5,10HNOT FOUND.) 

END 


SUBROUTINE MAXERC DATA, CHMAX) 
INTEGER DATA C 1 024) , CHMAX 
CHMAX =450 
DO 400 1=450,1024 
IFC DATA C I ) , GT.DATAC CHMAX) )CHMAX= I 
400 CONTINUE 

RETURN 
END 


SUBROUTINE SET1 CCHMAX, IDAT, POINT) 

C ROUTINE TO SET UP CALL TO FITALL TO FIND HEIGHT OF 

C NORMALIZATION PEAK 


INTEGER CHMAX, I DATC 1 024 ), DATADC 4096) , LIST C520) 

REAL POINTC100) ,PARAMC15) 

COMMON DATAD 

EQUIVALENCE C LI STC 1 ), DATADC 3073) ), CPARAMC I ), DATADC 4056) ) 

IMAX=CHMAX+20 

IMI N=CHMA X-20 

BG2 = IDAT C I MAX) 

BG1 = IDATCIMIN) 

DO 200 1 = 1 ,41 
I CO= IMI N+I - 1 
POINTCI)=IDATCICO) 

200 CONTINUE 

PARAMC2) = CBG1 ) 

PARA MCI ) = CBG2-BG1 ) /40 . 

PAP AMC5)=6./1 .665 
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PARAM(3) = FLOAT! I DAK CHMAX) )-PARAM(2) -PARAM(1 )*20. 
PAPA M (4) = FLOAT! CHMAX- 1 MIN+1 ) 

LIST(1)=41 
LI ST(2) = 5 
LI ST(3) = 1 5 
LIST(4)=0 
LIST(5)=5 
LIST(6)=5 
RETURN 
END 


SUBROUTINE OUTER! PARA M, LI ST, A NORM, ABNORM, CHMAX) 
INTEGER CHMAX, LI ST(520) 

REAL PARAMC1 5) 

PARAM(4)=PARAM(4)+FL0AT(CHMAX-21 ) 

IF(LIST(4) .NE.0) WRI TE (6 , 1 1 1 ) LIST(4) 

111 FORMAT! 1 X, 1 1 HERROR, TYPE ,12) 

ANORM=PARAM!3)*PARAM(5)*l .7725 
ABNOR M= 40 .*PARAM(2)+800.*PARAM(1 ) 

WRITEC6 , 123) P ARAM(4) , ANORM 
WRI TE (6 , 1 24) ABNORM 

124 FORMATC1 X, 1 5HBACKGROUND HAS ,F10.3,8 H COUNTS.) 

123 FORMAT!/, IX, 29HN0RMALI ZATION PEAK IN CHANNE 
1 L, F8 .2 ,4 H HAS, FI 0,2,8H COUNTS.) 

RETURN 

END 


SUBROUTI NE SUBTRC ANORM, APNORM, CHMAX, OUTPT, GAM) 

I NTEGER DATA (1 024) , CHMAX, DATADC4096) ,BKGND< 1 024) 

REAL PARAMC1 5) 

LOGICAL OUTPT, BITEST 
COMMON DA TAD 

EQUIVALENCE ( DATA! 1 ), DATADd 025) ), (BKGNDd ), DATAD(2049) ) , 
1 (PARAMO ), DA TAD! 4056) ) 

PARAM(4) =PARAM(4)+FLOAT(CHMAX-21 ) 
CN0RM=PARAM(3)*PARAM(5)*1 .7725 
CBNOR M= 40 .*PARAM(2)+R00.*PARAM( 1 ) 

WRITE(6, 123) PARAM(4) , CNORM 
WRI TE (6 , 1 24) CBNORM 

124 FORMATd X, 1 5HBACKGROUND HAS ,F10.3,8H COUNTS.) 

RAT=AB NORM/ CBNORM 
DO 5 1=3,1023 
IS=I-2 

DATA! I) = FLOAT! DATA! I > > -RAT* FL0AT(3 KGND! I- 1 )+BKGND( I ) + 

1 BKGND(I+1 )>/3. 

B KGND(I S) = FLOAT (DATA! I) ) +2.*R AT* FLOAT (B KGND! I -1 ) + 

1 BKGND! D+BKGND! 1+1 ) ) /3« 

IF(DATA(I) .LT.0) DATA(I)=0 
5 CONTINUE 

B KG ND(1 023) =0 
BKGNDd 024) =0 
CALL TOGGLS(IT) 

IF(BITEST(IT,15)) RETURN 
WRITE(6, 1 05) 

CALL I NPUT1 ( OUTPT, 5HOUTPT, 4) 

CALL INPUT1 (GAM, 5HGAMMA.2) 
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F0RMATC//,1 X, 2 9H NORMALIZATION PEAK IN CHANNE 
1L,F8.2,4H HAS, FI 0.2,8H COUNTS.) 

105 FORMATC1 X,25HD0 YOU WISH EXTRA OUTPUT?) 
RETURN 
END 


SUBROUTINE MULT( IMAX, I MI N, J) 

INTEGER CH0C30) , DA TAD (40 96) 

LOGICAL DOUBLE, UNDOUB, TRIPLE, UNTRIP, QUAD, UNQUAD, QUIN, UNQUIN 
COMMON/MUP/ DOUBLE, UNDOUB, TR IPLE, UNTRIP, QUAD, UNQUAD, QUI N , 

1 UNQUIN, JPEAK, 

1 GAM, I DUB, ITRIP, IQUAD, IQUIN 

COMMON DATAD 

EQUIVALENCE ( CH0C 1 ) ,DATAD(4021 ) ) 

DO UBLE=. FALSE. 

TRIPLE:. FALSE. 

QUAD=. FALSE. 

QUI N=. FALSE. 

I0:CH0C J) 

IMAX:CH0CJ)+I FIXC GAM) 

IMIN = CH0( J)-IFIX(GAM) 

IFC UNDOUB) GO TO 82 
IF( J.EQ. JPEAK) GO TO 82 

IFC CH0 ( J+l ) .LT.CI0 + IFIXC4.*GAM))) DOUBLE: . TRUE . 

IDUB=CH0( J+l ) 

IFC UNTRIP) GO TO 82 

IFC J.EQ. JPEAK-1 ) GO TO 82 

IF CCH0C J+2) .LT. C IDUB+IFIXC 4.*GAM) ) ) TRIPLE: .TRUE. 

TRIPLE: DOUBLE . AND . TR I PLE 
ITRIP:CH0C J+2) 

IFC UNQUAD) GO TO 82 

IFC J.EQ. JPEAK-2) GO TO 82 

IFC CH0C J+3) . LT.CITRIP+IFIXC4.*GAM))) QUAD = . TRUE. 

QUAD: TRIPLE. AND. QUAD 
IQUAD:CH0C J+3) 

IFC UNGUIN) GO TO 82 
IFCJ.EQ. JPEAK-3) GO TO 82 

IFCCH0C Jf4) ,LT.CIGUAD+IFIXC4.*GAM) )) QUIN:. TRUE. 

GUINsQUIN. AND. QUAD 
IQUIN=CH0C J+4) 

82 UNDOUB:. FALSE. 

UNTRIP:. FALSE. 

UNQUAD:. FALSE. 

UNQUIN:. FALSE. 

IFC DOUBLE) IMAX:IDUB+IFIXCGAM) 

IFC TRIPLE) IMAX = ITRIP+IFIXC GAM) 

IFCQUAD) IMAX=IQUAD+IFIXCGAM) 

IFCQUIN) IMAX = IQUIN+I FIXC GAM) 

RETURN 

END 


SUBRCUTI NE SI T2 C I MAX, I MI N , J , A NORM) 

C TO PREPARE FNIAL FIT TO DATA 

REAL PARAMC15),POINTC100) 

INTEGER DATA Cl 024) ,LISTC520) ,CH0 C30 ) , I RC 5) , DATADC 4096) 
LOGICAL DOUBLE, UNDOUB, TRIPLE, UNTR IP, QUAD, UNQUAD, QUI N , UNGUI N , 
1 BITEST 
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COMMON/MUP/DOUBLE, UNDOUB, TRIPLE, UNTRIP, QUAD, UNQUAD, QUI N, 
1 UNQUIN, JPEAK, GAM, IDUB, I TRIP, IQUAD, IQUIN 

COMMON DA TAD 

EQUIVALENCE C DATA( 1 ) , DATADC 1 025) ) , (LI STC 1 ) , DATADC 3073) ) , 

1 C CH0 ( 1 ) ,DATADC4021 ) ) , CPOINTC 1 ) , DATADC3593) ) , 

2 ( IRC I > , DATAD( 40 51 ) ) , CPARAMC 1 ) , DATADC4056) ) 
ITOT=IMAX-IMIN+l 

DO 84 1=1 ,ITOT 
IDOWN= IMIN+I- 1 
POINTCI)=FLOATCDATACIDOWN)) 

R4 CONTINUE 

I0 = CH0C J) 

PAR AM Cl 4) =FLOAT( IQUIN- IMIN+1 ) 

PARAMC 1 5) =GAM/1 .665 

PAR AM Cl 3) =FLOAT(DATA( IQUIN) ) 

PAP AM Cl I ) =FLOATC I QUAD- 1 MIN+1 ) 

PARAMC 1 2) = GAM/ 1 .665 

PARA MCI C) =FLOATCDATAC IQUAD) ) 

85 PARAMC7) =FLOATCDATAC ITRIP) ) 

PARAMC 8) =FLOATCITRIP-IMIN+l ) 

PARAMC 9) =GAM/1 .665 

86 PARAMC4) =FLOATCDATAC IDUB) ) 

PARAMC5) =FLOATC I DUB- IMIN+1 ) 

PARAMC6) =GAM/1 .665 

87 PAR AM Cl ) = FLOAT C DATA C 10) ) 

PARAMC2) =FLOATC I0-IMIN+1 ) 

PARAMC3) =GAM/1 .665 
LISTC1 ) = ITOT 

LI STC2) =3 
LISTC3) =15 
LI STC 4) =0 
LISTC5)=12 
LI STC6) =3 

IFC. NOT. DOUBLE) GO TO 88 
LI STC2) =6 
LI STC 6) =6 

IFC. NOT. TRIPLE) GO TO 88 
LISTC2) =9 
LI STC6) =9 

IFC. NOT. QUAD) GO TO 88 
LISTC2) =12 
LI STC6) =1 2 

IFC. NOT. QUIN) GO TO 88 
LISTC2) =1 5 
LI STC6) = 1 5 
LISTC5)=15 

88 CONTINUE 
NPARAM=LISTC6) 

CALL ZEROCIR) 

ISUM=0 

K=0 

DO 95 1 = 1 ,NPARAM,3 
K= K+l 

IFCPARAMCI) . LT.0.0000300*ANORM) IRCK)=1 
ISUM= ISUM+IRC K) 

95 CONTINUE 

LI STC2) = LI ST C2) -3* I SUM 
LI STC5) =LI STC5)-ISUM*3 
LI STC6) =LI ST C6) -3* I SUM 
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11=1 

DO 97 1 = 1, K 

IF( IR( I) . NE.0) GO TO 97 

96 PARAMdl )=PARAM(3*I-2) 

PARAMdl +1 ) = PA RAM( 3*1-1 ) 

PARAMdl +2) = PARA M( 3* I) 

98 11=11+3 

97 CONTINUE 

CALL TOGGLSC IT) 

IFC DOUBLE. AND. C. NOT. TRIPLE). AND. BITESTCIT, 14)) WRITE (6,1 15) 
IF( TRIPLE. AND. (.NOT. QUAD) . AND. BI TEST( I T, 1 4) ) WRITE(6,1 16) 
IF(QUAD.AND. (.NOT. QUIN) ,AND.BITEST( IT, 14)) WRITE (6,1 18) 
IF(QUIN.AND.BITEST(IT,14)) WRITE(6,119) 

RETURN 

115 FORMATO X,8H DOUBLET) 

116 FORMAT( 1 X,8H TRIPLET) 

118 FORMATO X, 1 1 H QUADRUPLET) 

119 FORMATO X, 1 1 H QUINTUPLET) 

END 


SUBROUTI NE IDENT( ELE, ENORM,B NORMX, DB, JPEAK, FNAME, NEXT) 

C SUBROUTINE TO MAKE FINAL IDENTIFICATION OF ELEMENTS, 

C AND MAKE CORRECTIONS FOR BETA PEAKS. 

REAL ELE(30) ,ENORM(30) ,B NORMX( 30 ) , DR(30) , FNAME(2) , 

1 ESI (24) , ES2(24) , AT OB (84) ,NAM(24) ,EFF(4) 

LOGICAL BITEST, BETA 

DATA ESI ( 1 ) , ES2 ( 1 ) , ATOB( 1 > , NAM( 1 )/2.62 ,2.81 ,0.,2HCL/, 

1 ESI (2) ,ES2(2) ,AT0B(2) ,N4M(2>/2.96,3. 19,0. ,2HAR/, 

2 ESI (3) ,ES2(3) ,AT0B(3) ,NAM(3)/3.31 3, 3. 589,0. 1 1 4 , 1 HK/ , 

3 ESI (4) ,ES2(4),ATOB(4),NAM(4)/3.691 ,4.012,0.0 90 9,2HCA/, 

4 ESI ( 5) , ES2( 5) , ATOB( 5) , NAM(5)/4.090 ,4 . 460 ,0 . 105,2HSC/, 

5 ESI (6) , ES2(6) ,AT0B(6) , NAM(6)/4. 5 1 0 ,4. 931 ,0. 105,2HTI/, 

6 ESI (7) ,ES2(7) ,ATOB(7) , N AM(7)/4 . 952 ,5 . 427 ,0 . 1 05 , 1HV/ 

DATA ESI (8) ,ES2(8) ,ATOB(8) , N AM(8 )/5 . 4 1 4 ,5 . 946 , 0 . 1 14,2HCP/, 

1 ESI ( 9) , ES2( 9) ,AT0B(9) , NAM(9)/5 . 898,6 .490,0 . 1 08, 2HMN/ , 

2 ESI (1 0) , ES2 (10) , ATOB( 10), NAM( 1 0 ) /6. 403 , 7.057 ,0. 1 06 ,2HFE/, 

3 ESI (1 1 ) ,ES2(1 1 ) , ATOB( 1 1 ) , NAM( 1 1 ) /6 . 930 ,7.649, 0 . 1 1 8 ,2HC0/, 

4 ESI (12) , ES2 ( 1 2) , ATOR( 1 2) , NAM( 1 2) /7 . 477 ,8.264 ,0 . 1 05 ,2HNI/ , 

5 ESI (1 3) ,ES2(13> ,ATOB( 13) , NAM( 1 3)/8.047 , 8. 904 , 0 . 1 1 3 ,2HCU/ 
DATA ESI (14) ,ES2(14),AT0B(14),NAM(14)/8.638,9.571 ,0 . 1 1 7 ,2HZN/ , 

1 ESI (15) ,ES2( 15) ,ATOB( 15) , NAMO 5) /1 0 . 543 , 1 1 . 725 ,0 .083,2HAS/ , 

2 ESI (16) ,ES2(16) ,ATOB( 16) , NAM(16)/1 1 .221 , 1 2. 495 , 0 . 1 3 ,2HSE/ , 

3 ESI (17), ES2( 17), ATOP (17) ,NAM(17)/1 1 . 923 , 1 3.290 ,0 . 1 36 ,2HBR/ , 

4 ESI (18) , ES2 (18) , ATOB( 1 8) , NAM( 1 8)/l 4 . 1 64 , 1 5 .834 ,0 . 1 51 ,2HSR/, 

5 ESI (19), ES2( 1 9) , ATOBO 9) ,NAM( 1 9) /I 5. 774, 1 7 .666 ,0 . 1 72 ,2HZR/ 
DATA ESI (20) ,ES2(20) , ATOB(20) , NAM(20) /I 7 . 478 , 1 9 . 607 ,0 . 1 8,2HM0/ , 

1 ESI (21 ) , ES2 (21 ) ,AT0B(21 ) ,NAM(21 )/8. 396 , 9 . 670 , 0 . 775 , 1 HW/ , 

2 ESI (22) , ES2 (22) ,ATOR(20) , NAMC22) /I 0.266, 12 .21 0 ,0 .672.2HTL/ , 

3 ESI (23) ,ES2(23) , ATOBC23) , NAM(23) /I 0 .549 , 1 2. 61 1 , 0 . 645 ,2HPD/ , 

4 ESI (24 ) , ES2 (24) ,AT0R(24) , NA M(24) /I 2. 966 , 1 6 .200 ,0 .553 ,2HTH/ 
DATA EPS, 1JNKNW/0 .060 ,2H--/ 

DATA EFF (1 ) , EFF(2) , EFF(3) , EFF(4)/1 .775,1.66,1 .66,2.00/ 

MEXT= NEXT- 1 

CALL BCDEXT( FNAME, MEXT) 

DO 4 1=3, JPEAK 
ELE ( I) =UNKNW 
DO 2 J=1 ,24 

IF( ARS( ENORMd) - ESI ( J) ) .LT. EPS) GO TO 1 
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GO TO 2 

1 BETA=. FALSE. 

DO 3 K= I, JPEAK 

IF(ARS(ENOPM(K) -ES2CJ)) ,GT. 0.090) GO TO 3 
IF(BNORMXCK) .GT.0.3*BNORMX(I)) RETA = . TRUE. 

B NOP MX ( K) =B NORMXC K) -ATOB( J)*BNORMX( I) 

IF( BNORMXC K) .LT.0.) BNORMXC K)=0. 

3 CONTINUE 

IF< J.LT .21 ) ELE ( I ) = NAM( J) 

IF( J.GE.21 .AND. BETA) GO TO 6 
GO TO 2 

6 ELE ( I) = NAM( J) 

BNORMX(I)=BNORMX(I)*EFFC J-20) 

DBC I) = DB( I)*EFF(H-20) 

2 CONTINUE 

4 CONTINUE 

CALL TOGGLSC IT) 

IFC , NOT.BITESTC IT, 1 3) ) GO TO 5 
CALL ENTER(3 , FNAME) 

WRITE (3,1 00) ( ELE ( I ) ,B NORMX( I), DP (I), 1 = 3, JPEA X) 

100 FORMAT ( 1 X, 1 A5 ,2 El 1 .5) 

CALL CL0SE(3) 

5 CALL TOGGLS(IT) 

IFC . NOT .BITEST( IT, 1 2) ) RETURN 

WRITE (6, 103) FNAME 
WRI TE (6 , 1 04) G 
FOR MAT ( 1 X,2HG= , F6 .4) 

WRI TE(6 , 1 01 ) 

WRITE (6, 102) (ELE (I) ,BNORMX(I) ,DB(I), 1 = 3, JPEAK) 
RETURN 

101 FORMAT ( 1 X, 7HELEMENT, 2X, 1 3HC0NCENTRATI0N, 9H ERROR) 

102 F0RMAT(4X,1A5,E12.6,E12.S) 

103 FORMATd HI ,5X,2A5) 

END 


SUBROUTINE ERRO(J,IMIN,DP) 

REAL ELE (30) ,B NORMX(30) ,DP(15),PARAM(15) 

INTEGER LI ST( 520) ,DATAD(4096) ,IR(5) ,CH0 (30) 

LOGICAL DOUBLE, UNDOUB, TRIPLE, UNTRIP, QUAD, UNQUAD, QUIN, UNQUIN, 
1 BITEST 

COMMON /MUP/DOUBLE, UNDOUB, TRIPLE, UNTRIP, QUAD, UNQUAD, QUIN, 

1 UNQUIN, JPEAK, GAM, I DUB, I TRIP, I QUAD 

COMMON DATAD 

COMMON/FI NOUT/ELE, A1 ,A2 

EQUIVALENCE ( LI ST(1 ) , DATA DC 3073) ) , (BNORMX( 1),DATAD(3861)), 

1 ( DB( 1 ) ,DATAD(3921 ) ) , ( IR( 1 ) ,DATAD(4051 ) ) , ( DATAD( 4021 ) , 

2 CH0 ( 1 ) ) , (PARAMO ) , DATAD( 4056) ) 

DATA RNONE/4 HNONE/ 

CALL TOGGLSC IT) 

IF(BI TEST (IT, 14)) WRITE (6, 111) LIST(4) 

111 FORMATd 1H ERROR, TYPE, 12) 

IF(. NOT. DOUBLE) GO TO 89 
IF( .NOT. TRIPLE) GO TO 90 
IF( .NOT. QUAD) GO TO 91 
IF(. NOT. QUIN) GO TO 92 
UNQUI N= .TRUE. 

RETURN 
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92 UNQUAD= .TRUE. 

RETURN 

91 UNTRIP=.TRUE. 

RETURN 

90 UNDOUB= .TRUE. 

RETURN 

R 9 PARAMCl ) =0 . 

PARAMC3) =0. 

PARAMC2) =FLOATC CH0C J)-IMIN+1 ) 
DP(1 ) = 0. 

DPC2)=0. 

DPC3)=0. 

RETURN 

END 


SUBROUTINE ENERGC IMI N, A1 ,A2,J,DP> 

REAL PARAMCl 5) ,ENORMC30> ,DPC15) 

INTEGER DATADC409S) , LI ST(520) , I R( 5) ,CH0 (30) 

COMMON DATAD 

EQUIVALENCE ( DATADC3R00) , ENORMC 1 ) ) , ( LI ST( 1 ) , DATAD( 3073 ) ) , 

1 (CH0C1 ) , DATAD C4021 ) ) , C IRC 1 ) , DAT AD ( 4051 ) > , 

2 CPARAMC1 ) ,DATADC4056)) 

K=J 

I SUM: 0 
DO 2 1=1,5 

IFCIR(I) .EG. 2) GO TO 2 
I SUM= I SUM+I R( I) 

2 CONTINUE 
K1 =0 

NPPEAK=LISTC6)/3 

NPEA K= NPPEA K+I SUM 

DO 3 1=1 , NPEA K 

MPEAK= NPEAK- 1+1 

M3PEAK=3* "PEAK 

IFC IRC MPEA K) . EG.0) GO TO 4 

JMPEAK=J+MPEA K-l 

PARA MC M3 PEAK) = 0. 

DPCM3PEAK) = 0 • 

PAPA MC M3 PEA K-l ) =FLOATC CH0 C JMPEA K) ) -FLOATC I MI N- 1 ) 

DPC M3PEAK-1 ) =0 ♦ 

PARA MC M3 PEA K-2) =0 . 

DPCM3PEA K-2) =0. 

GO TO 3 

4 N3PEA K=3* C NPPEA K-Kl ) 

PAPAMC M3PEAK) =PARAMC N3PEAK) 

DPCM3PEAK) =DPC N3PEAK) 

PAPAMC M3 PEA K-l ) = PARAMC N3PEAK- 1 ) 

DPC M3 PEA K-l ) =DPC N3PEA K-l) 

PAPAMC M3PEAK-2) =PAPAMC N3PEAK-2) 

DPC M3PEA K-2) = DPC N3PEAK-2) 

K1 =K1 +1 

3 CONTINUE 
LISTC6)=LISTC6)+3*I SUM 
NPARAM=LI STCS) 

DO 5 1 = 1 , NPARAM, 3 
IP2= 1+1 

PARAMCIP2)=PARAMCIP2)+FL0ATCIMIN-1 ) 

ENORMC K) =A1 *PARAMC IP2)+A2 
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K=K+1 

CONTINUE 

RETURN 

END 


SUBROUTINE LOCUM 

C TO LOCATE PEAKS IN X-RAY FLUORESCENT ENERGY SPECTRA.. 

C TO RE USED WITH IPK3, AS PART OF AN EXECUTE FILE. .SAVLOC 

INTEGER DATAC1 024) ,G4TE(5,30) ,CH0(30) ,DATAD(4095) ,CHMAX, 
1 CHI ,CH2,CH3,CH4,CH5,RKGND( 1024) , GAT1 (140) 

REAL C(36) ,B(36) ,PHI2,PHI 

LOGI CAL DOUBLE, UNDOUB, TR I PLE, UNTRI P, QUAD, UNQUAD, 

1 QUIN, UNPUI N 

COMMON/STUFF/GATE 

COMMON/MUP/DOUBLE, UNDOUB, TRIPLE, UNTRIP, QUAD, UNQUAD, 

1 QUIN, UNOUIN,JPEAK, GAM, IDUB, ITRIP, IQUAD,IQUIN 

COMMON DATAD 

COMMON/DIFF/C,R,M 

EQUIVALENCE ( DATA( 1 ) , DATAD (1 025) ) , ( CH0 C 1 ) , DATAD( 4021 ) ) , 

1 ( BKG NDC 1 ) , DA TAD (20 49) ) ,(GAT1 (1) , GATE (1,3)) 

CALL ZERO(GATl) 

NHMAX=470 

M=IFIX((0.6*GAM-1 .)/2.) 

IF( M. EG .0) M=1 

C COMPUTE COEFFICIENT CCJ) 

CALL ZERO ( C) 

CALL ZERO ( B) 

C(l)=-0.2 

C(2)=0.1 

C(3):0. 

DO 3 1 L= 1 , 5 
DO 30 J=1 ,36 
B( J)=C( J) 

DO 30 K= 1 , M 
K1 =J+K 
K2 - J— K 

IF(J-K.LE.P) K2 = K-J+2 
B( J) =B( J)+C( K1 )+C( K2 ) 

IF(K1 .GT.36) B( J) =0 

30 CONTINUE 

DO 31 J=1 ,36 
C( J) =B( J) 

31 CONTINUE 

C COMPUTE GENERALIZED SECOND DIFFERENCE, S 

C COMPUTE STANDARD DEVIATIONS IN S 

PHI=C( 1 )**2 
Ml = 5* Mf2 
DO 50 J=2 , Ml 
PHI=PHI+2.*C(J)**2 
50 CONTINUE 

PHI2= SQRT(PHI) 

C COMPUTE CRITERIA FOR PEAK LOCATION 

C NUEXPECTED NUMBER OF CHANNELS IN NEGATIVE PEAK 

Nl = IFIX(l ,22*(GAM+0 .5) ) 

N1 MAX=Nl+4 
N1 MIN=Nl-2 

C BEGIN TO SEARCH FOR PEAKS 

CH5=25 
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JPEAK=2 

58 JPEA K= JPEA K+l 

59 ICHNL=CH5 

60 IFCICHNL.EQ.NHMAX) GO TO 70 

ICHNL=ICHNL+1 

IFC SC I CH NL) .GE.-l . 6* PH I2*SRRTC FLOATC B KG NDC ICHNL-2)+l )) ) GO TO 60 
GO TO 62 

61 IF( ICHNL.EO. NHMAX) GO TO 70 

I CH NL = I CH NL+1 

62 IFC SC I CHNL+1 ) .LT.SC I CHNL) . OR. SC I CHNL- 1 ). LT. SC ICHNL) ) GO TO 61 
IF(SCICHNL) .LT.-l ,6*PHI2*SQRTC FLOATC BKGNDC ICHNL-2)+l ))) GO TO 63 
CH5= ICHNL 

GO TO 59 

625 C HA = ICHNL 

GO TO 64 

63 IFCICHNL.EQ.NHMAX) GO TO 70 

ICH NL= I CHNL+1 

64 IFCS(ICHNL) . LE.0 ,8*RC CH4) ) GO TO 63 
CH5 = ICHNL- 1 

ICHNL=CH4 

65 IFCICHNL.EP.l ) GO TO 59 

I CH NL= I CH NL- 1 

IFCSCICHNL) .LT.0.) GO TO 65 

CH3=ICHNL+1 

GO TO 67 

66 IFCICHNL.EO.l) GO TO 59 

ICHNL=ICHNL-1 

67 I F C SCI CHNL) . LT.PHI2* SQPTC FLO AT (P KGNDC ICHNL-2) +1 ) ) ) GO TO 66 
CH2= ICHNL 

GO TO 69 

68 IF(ICHNL.EQ.l) 30 TO 59 

I CHNL= I CHNL-1 

69 IFCSCICHNL) . GT.PHI2* SQRTC FLOATCB KGNDC ICHNL-2)+l )) ) GO TO 68 
CHI = I CHNL+1 

C COMPUTE EXPECTED WIDTHS 

N2 MIN = IFIXC CPHI2*SRRTC FLOATC DATAC CH4 ) +1 ) ) /ARSC SC CH4) ) ) 

1 *0.5* FLOATC N1 MIN) +0.5) 

N2 MAX = IFIXC CPHI2*SQRTC FLOATC DATAC CH4) +1 )) /A3SC SC CH4) ) ) 

1 *0.5* FLOA T C N1 MA X) +0 . 5 ) 

N3MIN = IFIXC FLOATC N1 MI N) * C 1 . -2 .* CPHI2*SQRTC FLOATC DATA C CH4) +1 ))/ 

1 ARSC SC CH 4)))) +0.5) 

N3MAX=IFIXC FLOATC N1 MAX)*C1 .-2.* CPHI2*S0RTC FLOATC DATA CCH4) + 1 ) )/ 

1 ABSCSCCH4))))+0.5) 


C IFCCH5-CH3+1 .LT.N1 MIN) GO TO 59 

C IFCCH5-CH3+1 .GT.N1 MAX) GO TO 59 

C IFC IE MAX.EO.0) GO TO 72 

C IFCCH3-CH2-1 .GT.N2 MAX) GO TO 59 

GO TO 73 

72 IFC CH3-CH2- 1 .GT. 1 ) CONTINUE 

73 IFCCH2-CH1 +1 .LT.N3MIN) CONTINUE 


CH0 C JPEAK) =CH4 
GATEC1 , JPEA K)=CH4 
GA TEC2 , JPEA K) = DA TA C CH4 ) 
GATEC3, JPEA K) =0 
GATEC4, JPEAK) = 128 
GATEC5 , JP EA K) =-10 
CALL SCAL1 C NX, NY, NF) 

GO TO 58 

70 JPEA K= JPEA K-l 

RETURN 
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TABLE I. - ELEMENTS RECOGNIZED BY SAVLOC 




(a) K series 


Element 

energy, keV 

Chlorine 

2.62 

Argon 

2.96 

Potassium 

3.313 

Calcium 

3.691 

Scandium 

4. 090 

Titanium 

4.510 

Vanadium 

4. 952 

Chromium 

5.414 

Manganese 

5.898 

Iron 

6.403 

Cobalt 

6.930 

Nickel 

7.477 

Copper 

8.047 

Zinc 

8.638 

Arsenic 

10. 543 

Selenium 

11.221 

Bromine 

11.923 

Strontium 

14. 164 

Zirconium 

15.774 

Molybdenum 

17.478 


TABLE II. - PARAMETERS USED TO 
DESCRIBE X-RAY EFFICIENCY 
FUNCTION F(E) 


I 

p(D 

1 

0.855 

2 

54. 874 

3 

-705.322 

4 

5499.011 

5 

-22895.21 

6 

48343.26 

7 

-40357.17 


(b) L series 


Element 

energy, keV 

L p energy, keV 

Tungsten 

8.396 

9.670 

Thallium 

10.266 

12.210 

Lead 

10.549 

12.611 

Thorium 

12.966 

16.200 
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RUN NUMBER SEPT 01 I 


NORMALIZATION PEAK IN CHANNEL 500.93 HAS 526R78.78 COUNTS. 
BACKGROUND HAS 152552.379 COUNTS. 

DNAMECA) = SEPT 
NEXT Cl) =0 


NORMALIZATION PEAK IN CHANNEL 501.48 HAS 648537.67 COUNTS. 
BACKGROUND HAS 56075.779 COUNTS. 

DO YOU WISH EXTRA OUTPUT? 

OUTPTCL) = F 
GAMMACR) = 4. 

ENERGIES FOR CALIBRATION: 

El CR) =6.403 

E2 CR) =12.611 

G=l .0000 


PEAK 

CHANNEL 

STRENGTH 

ENERGY 

CONCENTRATION 

CUGRAMS/IN**2) 

ERROR 

ELEM 

3 

70.035 

.000837 

3.701 

.29701+02 

.19374+01 

CA 

4 91.404 

DOUBLET 

.004320 

4.516 

.49B22+02 

.30408+01 

TI 

5 

128.004 

.000599 

5.912 

.22555+01 

.78361+00 

MN 

6 140.871 

QUA DR UP LET 

.077698 

6.403 

.221 15+03 

.24637+01 

FE 

7 

157.800 

.010517 

7.048 

,22002+02 

.15192+01 

FEBET 

8 

168.606 

.001 316 

7.461 

.23297+01 

.94050+00 

NI 

9 

183,969 

.004440 

8.047 

.64054+01 

.82467+00 

CU 

10 

1 99.256 

.032763 

8.630 

.3971 1+02 

.13505+01 

ZN 

11 

224.572 

.004366 

9.595 

.41648+01 

.61117+00 

ZNBET 

12 

249.423 

.029359 

10.543 

.23081+02 

.55729+00 

PBL 

13 

285.057 

.012955 

11.902 

.80868+01 

.28080+00 

BR 

14 

303.631 

.025692 

12.611 

. 1 4433+02 

.34989+00 

PBL 

15 321.212 

DOUBLET 

.001912 

13.282 

.97893+00 

.19976+00 

BRBET 

16 

359.355 

.003763 

14.737 

.16030+01 

.17434+00 


17 

368.690 

.001194 

15.093 

.48783+00 

.15608+00 


18 

385.632 

.003646 

15.739 

.13832+01 

.25890+00 

ZR 

1 9 

407.321 

.000938 

16.566 

.32513+00 

,71051-01 


20 

429.751 

.051486 

17.422 

.16323+02 

,29195+00 

MO 


Figure 5. - Sample of first page of teletype output. 


SEPT 

01 1 


ELEMENT 

CONCENTRATION 

ERROR 

CA 

.297007+02 

.193737+01 

TI 

.498224+02 

.304084+01 

MN 

,225552+01 

.783605+00 

FE 

.220910+03 

.246371+01 


.000000+00 

.151 91 9+01 

NI 

.232971+01 

.940501+00 

CU 

.640538+01 

.824675+00 

ZN 

.3971 09+02 

.135047+01 


. 000000+00 

.61 1 1 74+00 

PR 

.341605+02 

.824787+00 

BR 

.808675+01 

.280804+00 

-- 

.000000+00 

.349895+00 

- - 

.000000+00 

.199760+00 

-- 

.160305+01 

.174340+00 


.487833+00 

.156084+00 

ZR 

.138316+01 

.25 8902+00 


.325127+00 

.710506-01 

MO 

.163229+02 

.291 951+00 

Figure 6. 

- Sample of second page of 


teletype output. 
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RUN NUMBER SEPT 011 


NORMALIZATION PEAK IN CHANNEL 500.93 HAS 526878.78 COUNTS. 
BACKGROUND HAS 152552.379 COUNTS. 

DNAMECA) = SEPT 
NEXT (I) =0 


NORMALIZATION PEAK IN CHANNEL 501.48 HAS 648537.67 COUNTS. 
BACKGROUND HAS 56075.779 COUNTS. 

DO YOU WISH EXTRA OUTPUT? 

OUTPTCL) =T 
GAMMA (R) = 4. 

ENERGIES FOR CALIBRATION: 

El CR> =6.403 

E2 (R) =12.611 

G = 1 .0000 


PEAK CHANNEL STRENGTH 

ENERGY CONCENTRATION 

ERROR 



C UG RAMS/I N**2) 


3 70.035 

000837 

3.701 .29701+02 

.19374+01 

P< 1> = 213.043 

DELTA= 

6.546 


PC 2) = 70.035 

DELTA= 

.029 


PC 3) = 1.167 

DELTAS 

.038 


4 91.404 

004320 

4.516 .49822+02 

.30408+01 

PC 1)= 509.417 

DELTAS 

14.783 


PC 2) = 91 .404 

DELTAS 

.055 


PC 35= 2.521 

DELTAS 

.075 


DOUBLET 




5 128.004 

000599 

5.912 .22555+01 

.78361+00 

6 140.871 

077698 

6.403 .22115+03 

.24637+01 

PC 1>= 98.433 

DELTAS 

15.671 


PC 2) = 128.004 

DELTAS 

.235 


PC 3) = 1.807 

DELTAS 

.330 


PC 4 ) = 7328.194 

DELTAS 

44.861 


PC 5) = 140.87! 

DELTAS 

.012 


PC 6) = 3.152 

DELTAS 

.015 


QUADRUPLET 




7 157.800 

010517 

7.048 .22002+02 

.15192+01 

8 168.606 

001316 

7.461 ,23297+01 

.94050+00 

9 183.969 

004440 

8.047 .64054+01 

.82467+00 

10 199.256 

032763 

8.630 .39711+02 

.13505+01 

PC 1>= 959.873 

DELTAS 

32.626 


PC 2) = 157.800 

DELTAS 

.079 


PC 3) = 3.257 

DELTAS 

.110 


PC 4>= P 211.486 

DELTAS 

41.329 


PC 5>= 168.606 

DELTAS 

.280 


PC 6): 1 .849 

DELTAS 

.378 


PC 7) = 437.471 

DELTAS 

28.286 


PC B>= 183.969 

DELTAS 

.142 


PC 9) = 3.017 

DELTAS 

.188 


PC 1 0 > = 2772.324 

DELTAS 

48.293 


PCI1>= 199.256 

DELTAS 

.042 


PC 1 2) = 3.513 

DELTAS 

.057 


11 224.572 

004366 

9.595 .41648+01 

.61117+00 

PC 1)= 426.654 

DELTAS 

29.386 


PC 2) = 224.572 

DELTAS 

.157 


PC 3) = 3.042 

DELTAS 

.233 


12 249.423 

029359 

10.543 .23081+02 

.55729+00 

PC 1>= 2379. 97R 

DELTAS 

27.044 


PC 2>= 249.423 

DELTAS 

.029 


PC 3) = 3.667 

DELTAS 

.046 


13 285.057 

012955 

11.902 .80868+01 

.28080+00 

PC 1>= 1016.039 

DELTAS 

16.254 


PC 2>= 285.057 

DELTAS 

.042 


PC 3) = 3.790 

DELTAS 

,07 0 


14 303.631 

025692 

12.611 .14433+02 

.34989+00 

PC 1)= 1891.706 

DELTAS 

20.430 


PC 2>= 303.631 

DELTAS 

.031 


PC 3) = 4.037 

DELTAS 

.054 


15 321.212 

001 912 

13.282 .97893+00 

.19976+00 

PC 1>= 165.805 

DELTAS 

14.989 


PC 2) = 321 .212 

DELTAS 

.245 


PC 3) = 3.428 

DELTAS 

.385 


DOUBLET 




16 359.355 

003763 

14.737 .16030+01 

.17434+00 

17 368.690 

001 1 94 

15.093 .48783+00 

. 1 5608+00 

PC 1>= 294.605 DELTA= 

13.175 


PC 2)= 359.355 

DELTAS 

.150 


PC 3>= 3.797 

DELTAS 

.241 


PC 4> = 99.216 

DELTAS 

13.019 


PC 5>= 368.690 

DELTAS 

.429 


PC 6) = 3.579 

DELTAS 

.669 


18 385.632 

0036 46 

15.739 .13832+01 

.25890+00 

PC 1) = 242.082 

DELTAS 

17.146 


PC 2)= 385.632 

DELTAS 

.266 


PC 3) = 4.477 

DELTAS 

.516 


19 407.321 

000938 

16.566 .32513+00 

.71051-01 

PC 1) = 146.899 

' DELTAS 

14.928 


PC 2) = 407.321 

DELTAS 

.157 


PC 3)= 1 .898 

DELTAS 

.221 


20 429.751 

051486 

17.422 .16323+02 

.29195+00 

PC 1>= 3371.071 

DELTAS 

23.320 


PC 2) = 429.751 

DELTAS 

.025 


PC 3) = 4.540 

DELTAS 

.049 



Figure 7. - Sample of augmented page 1 teletype output. 


ELEM 

CA 


TI 


PIN 

FE 


FEBET 

NI 

CU 

ZN 


ZNBET 


PBL 


BR 


PBL 


BRBET 


ZR 


MO 
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Figure 8. - Flow 
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